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Abstract: We demonstrate how one can construct renormalizable perturbative ex- 
pansion in formally nonrenormalizable higher dimensional field theories. It is based 
on 1 /A'^^-expansion and results in a logarithmically divergent perturbation theory 
in arbitrary high space-time dimension. First, we consider a simple example of 
iV-component scalar filed theory and then extend this approach to Abelian and non- 
Abelian gauge theories with Nf fermions. In the latter case, due to self-interaction 
of non- Abelian fields the proposed recipe requires some modification which, however, 
does not change the main results. The resulting effective coupling is dimensionless 
and is running in accordance with the usual RG equations. The corresponding beta 
function is calculated in the leading order and is nonpolynomial in effective coupling. 
It exhibits either UV asymptotically free or IR free behaviour depending on the di- 
mension of space-time. The original dimensionful coupling plays a role of a mass 
and is also logarithmically renormalized. We analyze also the analytical properties 
of a resulting theory and demonstrate that in general it acquires several ghost states 
with negative and/or complex masses. In the former case, the ghost state can be re- 
moved by a proper choice of the coupling. As for the states with complex conjugated 
masses, their contribution to physical amplitudes cancels so that the theory appears 
to be unitary. 
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1. Introduction 

Nowadays it is popular to consider theories in extra dimensions as possible candidates 
for models of physics beyond the Standard Model |[l|, [^. Within the braneworld 
scenario one assumes that the matter fields are localized at the brane while the force 
carriers can travel in the bulk 0. Sometimes other fields might also live in extra 
dimensions. This means that one has higher dimensional QFT at least at short 
distances. However, it can hardly be considered as a consistent quantum theory 
beyond the tree level because of a lack of renormalizable perturbative expansion. 
Indeed, the usual coupling has a negative dimension, thus leading to power increasing 
divergencies which are out of control. In our previous work we studied the UV 
divergencies in scalar theories in extra dimensions within the perturbative expansion 
and demonstrated that although the leading divergences are governed by the one- 
loop diagrams even in the nonrenormalizable case, as was argued in this does not 
help to conquer them. 

Popular reasoning when dealing with extra dimensional theories relies on higher 
energy (string) theory which is supposed to cure all the UV problems while the low 
energy one is treated as an effective theory basically at the tree level. One way to 
do it is the Kaluza-Klein approach In this case, one takes the Fourier transform 
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over the extra dimensions and obtains an infinite tower of states witli quantized 
masses. Tlien one lias to sum over all these states. This sum is usually divergent 
and a special prescription is needed to regularize it. Doubtfully, however, that this 
approach solves the problem of nonrenormalizability in extra dimensional theories. 
As was shown in [0 , the properly renormalized four- dimensional theory never forgets 
its higher dimensional origin. It has an explicit cut-off dependence and can only be 
treated as an effective theory 

Here we make an attempt to construct renormalizable expansion in such formally 
nonrenormalizable theories using the well known technique of the expansion 



10], where in the scalar case N is the number of the scalar field components and in 



gauge theories it is the number of fermion flavours Nf. The number of colours Nc 
is kept fixed. This approach was successfully applied to non-linear sigma-models in 
3 dimensions |Tl| and to quantum gravity in four dimensions [|12| both of which are 
nonrenormalizable by power counting. Effectively, as we will show below, it leads 
to higher derivative theories and causes the usual problems of unitarity, locality 
and causality. However, these problems could be overcome though the analysis was 
performed mainly in the leading order JlT], |12]. 



We follow the approach of [11, 12] and apply it to theories in extra dimensions 
with the aim to construct renormalizable and unitary 1/A^ expansion suitable for 
perturbative calculations. We first consider scalar higher dimensional theories as an 
example [1^ and then treat the gauge theories with fermions in the same way [14 . 
The resulting perturbation theory is shown to be renormalizable, logarithmically 
divergent in any dimension D and obtains an effective dimensionless expansion pa- 
rameter. It is nonpolynomial in effective coupling, but polynomial in 1/A^ and obeys 
the usual properties of renormalizable theory. It might be either UV asymptotically 
free or IR free depending on the space-time dimension D. The original dimensionful 
coupling does not serve as an expansion parameter anymore and plays the role of 
mass which is also logarithmically divergent and multiplicatively renormalized. 

Within the dimensional regularization technique [jl5[ we performed the renormal- 
ization procedure in scalar and gauge theories in arbitrary odd space-time dimension 
and calculated a few terms of the 1 /N expansion. Even dimensions, in principle, can 
also be treated by this method; however, they lead to some complications due to the 
appearance of log terms. 

It is well known that the main problem of the 1/A^ expansion is to prove unitarity 
of a resulting theory since the analytical properties of the effective propagator change. 
When summing up the vacuum polarization diagrams to the denominator of a singlet 
field one gets an imaginary part and, in general, additional poles in the complex 
plane. These poles correspond to ghost states with the wrong metric and negative or 
complex masses. It is a common problem in any realization of the 1 /N expansion [|12|, 
|16|, . Effectively, it leads to higher derivative terms which may result in dynamical 
instability [0. This is another issue that we do not discuss here. Note, however. 
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the existence of "benign" quantum mechanical higher-derivative systems, where the 
classical vacuum is stable with respect to small perturbations and the problems 
appear only at the nonperturbative level |T^. The question of unitarity in higher 
derivative gravity in four dimensions was discussed in |1^, where the role of ghost 
states was emphasized. It was also shown that the higher derivative operators 
do not always improve the UV behaviour due to subtleties in analytical continuation 
from Minkowski to Euclidean metric. 

Below we consider the unitarity problem in detail and suggest a possible solution 
which seems to lead to a unitary theory in physical subspace. 



2. 1/N expansion. Scalar theory 

To illustrate the method, we start with the scalar theory. Let us take the usual N 
component scalar field theory in D dimensions, where D takes an arbitrary value 
(> 4), with the self-interaction. The Lagrangian looks like 

^ = ^M^-^^V^-^(0T, (2.1) 

where is the number of components of 0. We put into the normalization of 
the coupling so that A is fixed while N ^ oo. The theory is nonrenormalizable by 
power counting, the coupling A has negative dimension [\] = 2 — D/2. It is useful to 
rewrite eq. (|2.1|) introducing an auxiliary field a [0] 

C = lid,$r - \m^' - ^-(0^) + (2.2) 

Now one has two fields, one component and one singlet with triple interaction. 
Let us look at the propagator of the a field. At the tree level it is just "iA" , but then 
one has to take into account the corrections due to the loops of (j) (see Fig.l). 

* o --0--0-- 

Figure 1: The chain of diagrams giving a contribution to the a field propagator in 
the zeroth order of 1/A^ expansion 



If one follows the A^ dependence of the corresponding graphs, one finds out that 
it cancels: they are all of the zeroth order in 1/A^. Thus, one can sum them up and 
get ^ 

= - - - -^1-0---^ = l/A-n(p2)' (2.3) 
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where the polarization operator Il{p^) depends on D. In the massless case it looks 
like 

The obtained propagator ( |2.3|) has a typical for 1/iV expansion behaviour. Namely, 
it has a cut starting from = (for m = 0, otherwise from 4m^) and poles at nega- 
tive or complex depending on D. Notice that f{D) is finite for any odd D, despite 
naive power counting, and diverges for even D. This is due to the use of dimensional 
regularization: the one-loop diagrams in odd dimensions are finite since the gamma 
function has poles only at integer negative arguments and not at half-integer ones. 
This phenomenon can also be understood in using other regularization techniques. 
In general one has the UV divergence which has to be subtracted. This subtraction 
requires the redefinition of simple loop diagrams in D dimensions. However, the 
number of these diagrams is limited by [D/2] — 1, i.e., in 4 and 5 dimensions one 
has to define 1 diagram, in 6 and 7 dimensions - 2 diagrams, etc. In what follows 
we assume that this definition is made a'la dimensional regularization. Moreover, 
for simplicity of integration we limit ourselves to odd dimensions, which allows us to 
avoid the appearance of the log terms. 

A special issue is the existence of poles in the propagator. Usually they signal 
of the appearance of new asymptotic states which raise the problem of unitarity of 
a resulting theory. We address this problem in more detail in Sec. 7. Here we just 
mention that one can avoid poles on the real axis and have only complex conjugated 
pairs. This is enough for integration in Feynman diagrams. 

Thus, we have now the modified Feynman rules: the propagator is the usual 
one while the a propagator is given by eq. (|2.3| ). One can then construct the diagrams 
using these propagators and the triple vertex having in mind that any closed cycle 
of gives an additional factor of and any vertex gives 1 / \/N. 

Let us first analyse the degree of divergence. Let us start with the propagator. 
If the diagram with two external (p lines contains L loops, then it has 2L vertices, 
2L — 1 (f) lines and L a lines. Since each a line now behaves like l/p^~^, the degree 
of divergence is 

uj{G) = LD - {2L- 1)2 -L{D- A) = 2\ (2.5) 

for any D. Since this is a propagator, the divergence is proportional to and thus 
is reduced to the logarithmic one. 

Let us now take the triple vertex. If it has L loops, then one has 2L + 1 vertices, 
2L (f) lines and L a lines. Hence, the degree of divergence is 

u{G) = LD - {2L)2 - L{D - A) = 0\ (2.6) 

for any D. Thus, we again have only logarithmic divergence. 
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At last, consider the a propagator. In L loops it has 2L vertices, 2L (j) lines and 
L — 1 (J lines. The degree of divergence is 

uj{G) =LD- (2L)2 - (L - 1){D - 4) = D - 4. (2.7) 

This means that in odd D it has no global divergence (again we explore the properties 
of dimensional regularization) and the only possible divergencies are those of the 
subgraphs eliminated by renormalization of and the coupling. To see this, consider 
a genuine diagram for the cr-field propagator which is shown in Fig. 2, where the blobs 
denote the IPI vertex or propagator subgraphs. 




Figure 2: General type of the cx-field propagator 



After the R' operation^ we do not have any poles in the integrand for the remaining 
one-loop integral. What is left is the finite part containing logarithms of momenta. 
This final integration has the following form: 



ln"(A;V;u2) in™(A;2/p2) i^k^k^ik - p) 
k'^{k — pY 



where n, m, k are some numbers. We ignore here all the masses since they do not 
contribute to the UV behaviour. Due to the naive power counting of divergences in 
dimensional regularization we obtain the result proportional to T{2 — D /2) which is 
finite for any odd D. The logarithms can not change this property. 

To demonstrate how this works explicitly, we consider a particular example of 
the two-loop diagram. The result of the /^'-operation is shown in Fig. 3. 



R' 






Figure 3: Demonstration of the global divergence cancellation in the two-loop diagram 



^The R' operation means that we subtract from the diagram all divergent subgraphs 
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After subtracting the divergence in a subgraph we have prior to the last integration 

d'^-^'k T{D/2-l-e)V{2-e) 1 1 V{D/2-l) 

k^{p-ky^ ^ ^T{D/2-2)T{D/2 + l-2e){k^y eT{D/2-2)T{D/2 + iy 

The pole terms in the integrand cancel and expanding it over e one gets log(A;^). It 
is, however, easier to integrate it without expanding over e which gives 

T\D/2 -l-s)T{2- e)T{D/2 - 1 - 2e)T{2 + 2e- D/2) {p^l^-^ 
^ V{D/2-2)V{D/2 + l-2e)T{l + e)V{D-2-?,e) {p^f' 

1 T{D/2 - 1) T\D/2 - 1 - e)T{2 + e - D/2) ^ 

e T{D/2-2)T{D/2 + l) T{D -2-2e) {p^Y ^ ^' 

Thus, after the R' operation the diagram is finite and we do not need the a field 
renormalizat ion . 

This way one gets the perturbative expansion with only logarithmic divergences. 
This is not expansion over dimensionful coupling A but rather 1/N expansion with 
dimensionless parameter. 



3. Properties of the 1/N expansion 

Consider now the leading order calculations. We start with the terms for the 
propagator of (f) and the triple vertex. One has the diagrams shown in Fig.4. Notice 
that besides the one-loop diagrams in the same order of the 1/N expansion one has 
the two-loop diagram for the vertex. 



Figure 4: The leading order diagrams giving a contribution to the (f) field propagator 
and the triple vertex in expansion 



Let us start with the diagram a). One has 

J {2n)DN [{k - pY - m2] [1 /A - Yi{k'^)] ' ~ ~ 

Since we are interested in the UV behaviour we can omit the mass from the field 
propagator and "1/ A" from the a field propagator and take the massless limit of the 
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polarization operator 11 (A;^) . Wc will restore them when discussing the analytical 
properties. Then the UV asymptotics is given by 



a 



{27r)^Nf{D) (A; -p)2(-A;2)^/2- 



One can see that the original coupling A plays the role of inverse mass and drops out 
from the UV expression. What is left is a dimensionless 1/N term. 

Calculating the singular parts of the diagrams of Fig.4 in dimensional regular- 
ization with D' — D — 2e one finds 

Diaq.a ^ —rzA, Diaq.h =^ —r;B, Diaq.c ^ —r;C, (3.1) 
A = 2r(D - 2) B- ^ A C- ^^^-^\ \ 

r(D/2-2)r(D/2- i)r(D/2 + i)r(2-D/2)' 4-l> ' a-d ' 



The corresponding renormalization constants in the MS scheme then are 

Z.-' = 1 - j^, (3.2) 

There is no any coupling in these formulas, its role is played hj 1/N which is therefore 
infinitely renormalized. This seems to be unsatisfactory and to overcome this problem 
we introduce a new dimensionless coupling h associated with the triple vertex (and 
not with the a propagator) as 

r ^ I2 



Then in the leading order in 1/A^ the renormalization constants and the coupling 
take the form 

Z-i = 1 - (3.4) 
7 , hB H'C 

(3.6) 

This is not, however, the final expression. To sec this, we consider the next 
order of the 1/N expansion. The corresponding diagrams for the (j) propagator are 
shown in Fig. 5. Again one can see that the l/N"^ terms contain not only the two- loop 
diagrams but also the three- and even four-loop ones. 
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Figure 5: The second order diagrams giving a contribution to the field propagator 
in expansion 



All these diagrams are double logarithmically divergent, i.e., contain both single 
and double poles in dimensional regularization. We calculate the leading double pole 
after subtraction of the divergent subgraphs, i.e. perform the i?'-operation. The 
answer is: 

-^^ACh^ Diag.e -^^A^ Diag.f --l^^ABh', 
D^^9-9 ^ -^ACh\ (3.7) 

Here we face a problem, namely, subtracting the divergent subgraphs in the graphs 
e-g, we get the diagram which is absent in our expansion, since it is already included 
in our bold a line (see Fig. 6). 

-O- 

/ \ 

f % 
I t 

I I 

Figure 6: The "forbidden" loop diagram 
There would be no problem unless this diagram is needed to match the so-called pole 



equations |^ which allow one to calculate the higher order poles in the Z factors from 
the single one. However, if we include this diagram in the a line, it will not change the 
latter, except for the additional h factor coming from the vertex and not compensated 
by the propagator. Apparently, one can continue this insertion procedure and add 
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any number of such loops not changing the order of expansion. The result is 
the sum of a geometrical progression 



which should multiply every a line. Altogether this leads to the following effective 
Lagrangian for UV perturbation theory 

= \{d,$f - + + \f{D)a{dr''-Ml + h). (3.8) 



Having all this in mind we come to the final expressions for the Z factors within 
the 1/N expansion: 

I f Bh C/i^ \ 1 
^1 I Ah 1/3 A^h'^ ABh'^ 2 A^h^ 



^ eNl + h e^N^ \2 (1 + hf [l + hf ?>{! + hf 

2 ABh^ 4 ACh"^ ACh^ \ 1 x , . 



4. expansion. QED 

Let us consider now the usual QED with Nf fermion fields in D dimensions, where 
D takes an arbitrary odd value. The Lagrangian looks like 

1 1 _ . _ e - - 

£ = —A^^,A^ - d^A^f - —{d^A^f + iV'.^V'i - m^i%lji + ^=ilj,AiPi. (4.1) 
4 za V^/ 

According to the general strategy, we now have to consider the photon propa- 
gator. Since due to the gauge invariance the polarization operator is transverse, it 
is useful to consider a transverse (Landau) gauge. This is not necessary but simpli- 
fies the calculations. Then in the leading order of the 1/N expansion one has the 
following sequence of bubbles (see Fig. 7) 




Figure 7: The chain of diagrams giving a contribution to the A field propagator in 
the zeroth order of the 1/A^/ expansion 
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summed up into a geometrical progression. This is nothing more than the renormalon 
chain |22|. The resulting photon propagator takes the form similar to that for an 
auxiliary field o in scalar case 

DM = (a"" - ^) i + e2jp)(_p2)D/2-2' (4-2) 



where 



^^^^ r2(D/2)r(2-D/2) 



2D-[D/2]-lrp)^D/2 

and we put m = for simplicity. 

This practically coincides with the expression obtained in scalar theory and all 
the following steps just repeat those in the latter. We change the normalization of 
the gauge field — > A^/e and introduce the dimensionless coupling h associated 
with the triple vertex, so the effective Lagrangian takes the form 



+iilJidilJi - mi/ji'ilJi + -^^^^ipiAil)i. (4.3) 

This new dimensionless coupling h enters into the gauge transformation and plays 
the role of a gauge charge. The old coupling e, on the contrary, is dimensionful and 
acts as a mass parameter in a gauge propagator. Since the coupling constant h is 
dimensionless the effective Lagrangian (|4.3|) when omitting the first term is conformal 



as considered in p3| where the theory was taken in D = 3. 

Again, one has the modified Feynman rules with the photon propagator that 
decreases in the Euclidean region like l/(p^)^/^~^, thus improving the UV behaviour 
in a theory. The only divergent graphs are those of the fermion propagator and the 
triple vertex. They are both logarithmically divergent for any odd D. The photon 
propagator is genuinely finite and may contain divergencies only in subgraphs. One 
basically has the same graphs as in a scalar theory but with solid lines being the 
fermion ones and the dashed lines being the photon one. 

The only difference (or simplification) comes from the Furry theorem and the 
gauge invariance. Namely, all triangles with three photon external lines vanish due to 
the Furry theorem and the gauge invariance which connects the fermion propagator 
with the triple vertex implies that Zi = Z2. This relation holds in the 1/iVj expansion 
hke in the usual PT. Thus, using the notation of a previous section, in the leading 
order one has 

r(D)(D-i)(2-D/2) 

2\D/2]+iY{2-D/2)V{D/2 + l)T\D/2y ' 



The same results were obtained in |^ where the author calculated the anomalous 
dimensions at the D-dimensional critical point where the fields obey asymptotic 
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scaling and are conformal. This leads to the following renormalization constants in 
the leading order in l/Nf-. 

and, consequently, hs = h. Hence, in odd-dimensional QED in the leading order of 
the 1/iV/ expansion one does not need the coupling constant renormalization; only 
the wave function renormalization remains. This means that the coupling is not 
running. 

In the second order one again has the same diagrams as in a scalar theory but 
with vanishing triangles. The renormalization constant in the second order is also 
essentially simplified compared to the scalar case and looks like 

I Ah 11 ^2/^2 1 , , , 

Like in the scalar case the original dimensionful coupling e is not an expansion 
parameter anymore, but plays a role of a mass and is multiplicatively logarithmically 
renormalized. The leading order diagrams are shown in Fig. 8. 

■~^^y— ™0™ 

a b c 

Figure 8: The first order diagrams giving a contribution to the 1/e^ renormalization 
in the 1/iVj expansion 

They give the following contribution: 

D^ag.a ^ ^^^^^^,F, D^ag.h ^ ^^^^^^,E, Drag.c 0, (4.7) 

r(D + l)iD/2 - 1){D - 1)2(2 - D/2) D^ + D/2-9 

2i?/2+ir(2-D/2)r(D/2 + 2)r2(D/2)' d/2{d/2-i){d -i) • 

So one has 

5. 1/N/ expansion. QCD 

Consider now a non-Abelian theory with Nf fermions. Notice that in QCD, contrary 
to QED, all Feynman diagrams contain group factors so that the actual expansion 
parameter becomes N^/Nf, thus requiring that this ratio is small. At the same time. 
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to preserve asymptotic freedom in 4 dimensions one needs Nc/Nf > 2/11. So one has 
some interval where the Nc/Nf expansion might be vahd. Of course, in non-Abehan 
theories the 1/Nc expansion would be preferable, since it accumulates the interactions 
of the gauge fields, however, in this case already the lowest approximation consists 
of all planar diagrams and is not known . 

In the non-Abelian case one has some novel features due to the presence of the 
triple and quartic gauge vertices and the ghost fields. Similar to ( [l.lD we write down 
the Lagrangian for the gauge fields and Nf fermions as 



where 



pa 
/if 



d,Al - d^Al + 



IN 



■■f 



f 



abc Ab AC 



D, 



A" 



Like in QED we choose the Landau gauge and sum up the fermion bubble diagrams 
into the denominator of the gauge field propagator 



ab 



i6° 



{9 



p2 l+g^f{D){-p^)^/^-^' 



(5.1) 



where the coefficient f{D) differs from the Abelian case only by the colour factor 
T(i?) 

r2(D/2)r(2-L)/2), 



;T{R) 



2D-[d/2]-iy(^D)t^d/2- 

and again we put m = for simplicity. 

In the non-Abelian case, contrary to the Abelian one, one has the triple and 
quartic self-interaction of the gauge fields. These vertices, which are suppressed by 
and 'i^/Nf, respectively, obtain loop corrections of the same order in 1/Nf. 
The effective vertices in the leading order are given by the diagrams shown in Fig. 9 
and 10. 





Figure 9: The diagrams giving a contribution to the A^ term in the zeroth order 
of the 1 / Nf expansion 
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Figure 10: The diagrams giving a contribution to the term in the zeroth order 
of the 1/A^/ expansion 

Thus, besides the modification of the gauge propagator one has the modified 
vertices. The effective Lagrangian in the case of vertices is not given by a simple 
local expression due to complexity of the loop diagrams. So we keep it in the form 
of the diagrams which have to be evaluated in integer dimension. Due to the rules 
of dimensional regularization they are finite for any odd D, otherwise one has to 
redefine them. What is crucial, however, is that there are only three diagrams which 
have to be redefined. Hence, after rescaling the gauge field ^fi/ 9 one obtains 

the following effective Lagrangian: 

^^// = - +xix + XX ) (5-2) 



^otg \/Nf 



Notice that dimensionful coupling g drops from all terms except for the first one and 
is not an expansion parameter anymore. 

Calculating the degree of divergence after summing up the diagrams of the zeroth 
order, similar to the scalar case and QED, one has only four types of logarithmically 
divergent diagrams: the fermion and the ghost propagators, the fermion-gauge- vertex 
and ghost-gauge-ghost vertex. The gauge propagator as well as pure gauge vertices 
are finite and may contain only divergent subgraphs. 

The next step is the introduction of a dimensionless coupling h. Here one should 
be accurate since this coupling enters not only into the triple gauge-fermion vertex, 
but due to the gauge invariance should be present in gauge and gauge-ghost vertices. 
It should be the same in all three of them. In the case of a gauge theory, the coupling 
h enters the gauge transformation and acts as a gauge charge of the fermion and gauge 
fields. 

When constructing the Feynman diagrams, one reproduces the one-loop cycles 
that are already present in the effective Lagrangian ( |5.2| ) but with additional factors 
h. In the scalar or QED case, this happened only for the propagator, but here it is 
also true for the vertices. As a result, the final expression for the effective Lagrangian 
takes the form 



^ff = -—.^PP' - (HJ- + ^U^ ) (1+^) (5-3) 
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where 



\/h \/h 



Consider now the leading order calculations. We start with the 1/Nf terms for 
the fermion and the triple fermion-gauge-fermion vertex. The diagrams are shown 
in Fig. 11. The first two are the same as in QED. The third diagram contains new 
effective vertex which includes the usual triple vertex and the fermion triangle. The 
usual vertex does not give a contribution since it is finite by a simple power counting. 
At the same time, the fermion triangle is momentum dependent and the resulting 
diagram is logarithmically divergent. 





Figure 11: The leading order diagrams giving a contribution to the ip field propagator 
and the triple vertex in 1 / Nf expansion 

Calculating the singular parts of the diagrams of Fig. 11 in dimensional regular- 
ization with D' = D — 2e one finds 

Diag.a A, Diag.h B, Diag.c C, 

eNf 1 + h eNf 1 + h eNj 1 + h 

T{D){D-1){2-D/2)Cf 

2[D/2]+ir(2-D/2)r(Z}/2 + i)r2(D/2)T' ^ ' 

- Ca/2 , il-D/2)CA , 

^~ C; ^~~2{2-D/2)Cf' 



which is again in agreement with Notice that the third diagram is proportional 
to h/{l + h) instead of h"^/ (1 + h)'^ as in the scalar case. The reason is that now we 
have an effective triple gauge vertex proportional to Vh{l + h) instead of h that 
cancels one factor oi h/{l + h). 

Therefore, in the leading order in the 1/iV/ expansion the renormalization con- 
stants take the form 

^2"' = 1 - ^T^^. (5.5) 
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To check the gauge invariance, we calculated the renormalization of the coupling 
through the gauge-ghost interaction. The leading diagrams are shown in Fig. 12. 




Figure 12: The leading order diagrams giving a contribution to the ghost field propagator 
and the triple vertex in 1 / Nf expansion 

Calculating the singular parts of the diagrams in dimensional regularization one finds 

1 h . 1 h 1 h 

Diag.a ^ —-——A , Diag.h =^ B, Diag.c C, 

eNf 1 + h eNf 1 + h eNj I + h 



2[^/2l+2r(2 - D/2)r(Z}/2 + l)r2(D/2)T' 
which gives the following renormalization constants in the ghost sector 



Zi = 1, (5.9) 
~ T 1 A'h , , 

= ' - 7n;ttt,- (^■^°' 

One can see that the following relation holds: 

A + 5 + C = A' + S' + C", (5.12) 

which follows from the gauge invariance. 

We look now at the next-to-leading order to compare it with the scalar case. 
The corresponding diagrams for the fermion propagator are shown in Fig. 13. They 
require some explanation. The first line of diagrams in Fig. 13 is obtained from the 
one-loop diagrams of Fig. 11 by inserting into the vertex or the fermion line of the one- 
loop divergent subgraphs from Fig. 11. For example, the diagram d in Fig. 13 is the 
diagram a from Fig. 11 with divergent one- loop subgraph c from Fig. 11 substituted 
instead of the initial vertex. The second line of the diagrams in Fig. 13 is obtained 
from the "forbidden" diagram of Fig. 14 by inserting the same one-loop divergent 
subgraphs from Fig. 11. The diagram e is the diagram of Fig. 14 with insertion of the 
subgraph a from Fig. 11 into the fermion line (see Fig. 15) and the diagram g comes 
from the insertion of the subgraph c from Fig. 11 instead of one of the vertices in the 
fermion loop (see Fig. 16). 
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A 



a 





f 



Figure 13: The second order diagrams giving a contribution to the fermion propagator 
in the 1/Nf expansion 




Figure 14: The "forbidden" loop diagram 




Figure 15: The diagram e from Fig. 13 as a result of insertion of the diagram a from 
Fig. 11 into the fermion line. 




Figure 16: The diagram g from Fig. 13 as a result of insertion of the diagram c from 
Fig. 11 instead of one of the vertices in the diagram from Fig. 14. 

All these diagrams are double logarithmically divergent, i.e., contain both single 
and double poles in dimensional regularization. We calculated the leading double 
poles after subtraction of the divergent subgraphs, i.e., performed the i?'-operation. 
The answer is: 

1 A^h^ , 1 ABh^ 1 A^h^ 



Diag.a 



e^N] 2(1 + hf 



Diag.h 



e^Nj (1 



Diag.c 
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Diag.d 
Diag.f 



1 2ABh^ 



Diag.e 



Diag.g 



(5.13) 



2ACh^ 



6^Nj3{l + hf' ''''' emj3{l + hY' 

We performed also the calculation for the fermion-gauge-fermion vertex but do not 
present the diagram-by-diagram result because of the lack of space and give only the 
final answer. 

Having all this in mind we come to the final expressions for the Z factors in the 
second order of the 1/iVj expansion in the fermion sector: 



1 - 



1 {B + C)h 



+ 



eNf 1 + h 

2{B + Cfh^ 2A{B + C)h^ 

3 + 3 

1 Ah 1 



3{B + Cfh^ A{B + C)h^ 

2 {1 + hy + {i + hy 



1 - 



+ 



3 A^h^ 



0{ 



1 



(5.14) 



eNf 1 + h 
2A{A + B + C)h^ 
'3~ 



e^Nj y2{i + hy 



+ 



A{B + C)h^ 

{i + hy 



+ 0{ 



(1 + hy 

The same calculation in the ghost sector gives 
Zi = 1, 



(5.15) 



(5.16) 



1 A'h 

~ sNfl + h^ 



em] 



3 A'^h^ 
2(l + /i)2 



2A'{A + B + C)h^' 
3~ 



Of 



Notice the absence of the ghost-gauge-ghost vertex renormalization. 

The final second order expression for the coupling renormalization calculated in 
both ways having in mind relation (|5.12|) is 



1 2{A + B + C)h 



eN 



f 



1 + h 



(A + B + Cfh^ 
{1 + W 



A {A 
3 



^B + CfK^ 
(1 + hf 



+0{ 



(5.17) 

Like in the scalar and QED case, one can also calculate the renormalization of 
the original coupling g"^ . The leading order diagrams are shown in Fig. 17 which 



eNf 



give the following singular parts like in [25 




a bed 

Figure 17: The first order diagrams giving a contribution to the l/g"^ renormalization 
in 1/A^/ expansion 
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^-^-^i^(r|^^' ^-^-^ ^ i^(r^^' (5-18) 

r(D + l){D/2 - 1){D - 1)2(2 - D/2)Cf 



F 
E 
G 



2P/2]+2r(2 - D/2)V{D/2 + 2)V^{D/2)T{R) ' 

D/2(D/2- 1)(Z}-1) Cf 
A{D/2f - Q{D/2f + 18(D/2)4 - Q7{D/2)^ + 85{D/2f - 19D + 6 Ca 
2{D - 1)2(1 - D/2)2(2 - D/2)D C^^ 



D^-D^/2-2D + l Ca 
D{1- D/2){2- D/2){D~iyCF ' 
The corresponding renormalization constant looks like 

7 1 ^ {F + E + G + H)h' 

6. Renormalization group in 1/N expansion 

Having these expressions for the Z factors one can construct the coupling constant 
renormalization and the corresponding RG functions. One has as usual in the di- 
mensional regularization 

hs = if^rhz^z^' = {f^y (h + Y, ^11^] , (6.1) 



n=l 



n=l 

where the first coefficients a„ and can be deduced from the Z factors. 

This allows one to get the anomalous dimensions and the beta function defined 

as 

7(/i, N) = -/i2^ log Z = /i J-ci, (6.3) 

/3(/i,iV) = 2/1(71 + 72) = (/i-^-l)ai. (6.4) 

ah 



We first consider the scalar case. With the help of eqs. (p.9| , |3.1CI| ) one gets in the 
leading order of 1/A^ expansion2 

1 Ah , I f Bh 2Ch^ \ , , 



^Notc that the anomalous dimension of a field 72, is defined with respect to Z2 ^. 
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n,. .r^ 1 f2(A + BW ACh^ , , , 



It is instructive to check the so-called pole equations |^ that express the coef- 
ficients of the higher order poles in e of the Z factors via the coefficients of a simple 
pole. For ^2"^ one has, according to (|3.1CI|) , 

1 Ah 

C2{h, N) = — { -7-^ + 7-^ + 



m \2{i + hy (i + /i)2 3{i + hy 

2 ABh^ 4 ACh^ ACh^ \ 
^3(l + /i)3 + 3(l + /i)3 + (l + /i)V ■ 

At the same time the coefficient C2 can be expressed through ci via the pole equations 
as 

which gives 

^dc2 _ 1 Ah Ah ^ I f2{A + B)h'^ , ACh^ \ A 



dh m{i + hyi + h \ (1 + hy (i + /i)^ ; (i + h) 



Integrating this equation one gets for C2 the expression coinciding with ( |6.8|) which 
was obtained by direct diagram evaluation. Notice that to get this coincidence the 
/i-dependence in the denominator of eqs.( ^.9| , |3.10D was absolutely crucial. 

We have also checked the pole equations for the renormalized coupling. They 
look as follows 

(AA-l)a„^,%l. (6.10) 

In the leading order in h when 



one should have a geometric progression 

a„(/i,iV) = ai(/i,iV)". 

We have checked this relation up to three loops and confirmed its validity. 

Having expression for the (3 function one may wonder how the coupling is run- 
ning. The crucial point here is the sign of the f3 function. One has 

, / 2/^2 D(D - 3)h^ 
AT{D - 2) 7—^ + ^ ^ 



(3(hN) = -- + (^ + ^)^ / (6 11) 

' ^ Arr(D/2-2)r(D/2-i)r(D/2 + i)r(3-L'/2)' ^ ' 
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It can also be rewritten as 

that clearly indicates that the theory is UV asymptotically free for D = 2k + 1, k - 
even and IR free for fc-odd. Solution of the RG equation looks somewhat complicated, 
but for the small coupling in the leading order it simply equals the usual leading log 
approximation 

For example, for D = 5,7 the beta function equals (3o = — 256/157r^A^ and 2^^/1057r^A^, 
respectively. 

We now come to the gauge theories. With the help of eqs. (|5.14| , |5.15|) one gets in 
the leading order of the expansion 

n .r . 1 , 1 {B + C)h 

72(ft,iV,) = -^^^, ^,(h,Nf)^0(j^), (6.15) 

The situation is similar to that in scalar theory. Only the value of coefficients 
are different. This, however, does not influence the pole equations. They remain to 
be vahd. 

Equation ( p. 17] ) gives us the sign of the beta function. In the leading order one 

has 

^^.(h) = rp)p-i)c^ ^ 

dt ' 2[^/2]+2r(2-D/2)r(D/2 + i)r2(D/2)Ar^r(i + /i)2' ^ ^ 

which means that contrary to the scalar case ( |6.12|) (3{h) > for D = 5, j3{h) < 
for D = 7 and then alternates with D as in the scalar case. 

Solution to eq.( |6.17| ) for small h is again reduced to the usual one. As for the 
original couplings, there is no simple solution either except for the QED case, where 
the coupling h is not running and solution of the RG equation for 1/e^ with fixed h 
is 

1-1 (^^ ^ 

e2 el \pI 

with the anomalous dimension 



7 



(6.18) 

T{D){D - l){D/2 - 2){D - 3){D + 2){D - 6) 

2[i?/2]+2rp/2 + 2)r2(D/2)r(2 - D/2)Nf (1 + /i)3 ■ 
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The sign of 7 depends on D. For D = 5, 7 7 > 0, for D = 9 7 < and then 
alternates with every odd D. Eq.( |6.18D reminds the power law behaviour of the 
initial coupling in extra dimensions within the Kaluza-Klein approach though 
anomalous dimension 7 is different. 



7. Analytical properties and unitarity 

Consider now the analytical properties of the propagator and related problem of 
unitarity. The problem is common to scalar and gauge theories so for simplicity we 
concentrate on the sigma field propagator (|2.^). Besides the cut starting from 4m^ 



it has poles in the complex p plane. Hence, knowing the analytical structure, one 



can write down the Kallen-Lehmann representation ||28 
Let us first consider the massless case 



Depending on a sign of f{D) there are two possibilities: either one has a pole at real 
axis and (possibly) pairs of complex conjugated poles (/(-D) < 0, D=5,9,...) or one 
has only pairs of complex conjugated poles (/(-D) > 0, D=7,ll,...) and all the rest 
appears at the second Riemann sheet. We consider the cases of D = 5 and D = 7 as 
the nearest options. One has, respectively, 

D^iP) = -Zr^77^^Tr-7TT^ + - / ...... , (7-2) 



and 



Drip' 



p2 + (2567r/A)2 TT Jq p2 — (2567r)2 + A^m^ ' 



^ 81927r^ ■j2/3^7rr/3 2 81927if,^2/3g-7ri/3 



3 p2 _^ (8]i^)2/3g-27ri/3 3 p"^ + SI^tH ^2/3g27ri/3 

dm^ 81927r2A2rm2)3/2 



n Jo p2 _ ^2 A2(r7i2)3 + (81927r2)2 ' 



(7.3) 



Notice that the continuous spectrum has a positive spectral density and corresponds 
to production of real pairs of fields (or pairs of fermions in the gauge case). These 
states are present in the original spectrum and cause no problem with unitarity. This 
analysis was performed at the tree level in |]TT| and can be extended to any number of 
loops. One can show that all the cuts imposed on diagrams when applying Cutkosky 
rules m in any order of perturbation theory lead to the usual asymptotic states on 



mass shell and no new states appear. 

The problem comes with the poles. One can see that the pole terms come with 
negative sign and, therefore, correspond to the ghost states [|^. For D = 5 one has 



only one pole at the positive real semiaxis while for D = 7 one has a pair of complex 
conjugated poles, as shown in Fig. 18. 
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Figure 18: The analytical structure of the auxiliary field propagator 



The presence of these ghost states is the drawback of a theory. They signal 
of instability of the vacuum state. Indeed, as it was shown in [^], the vacuum 
might be unstable with respect to appearance of condensates. This will lead to 
additional diagrams similar to those in QCD. However, they do not seem to improve 
the situation. Thus, one has either to try to get rid of ghost poles or to make sure 
that they do not give a contribution to physical amplitudes. 



Let us first see what happens if one takes a nonzero mass of the 
polarization operator then is 



field. The 



T{2-D/2) 

2D+lj^D/2 



dx{—p'^x{l — x) + m 



2\D/2-2 



(7.4) 



For D = 5 one has 



-m 



4V^+(a-4)ln(f^^ 



(7.5) 



2 

where a = Since the existence of a pole is governed by the equation 



n(/) = i/A, 



one has to check whether this equation is satisfied somewhere in the complex 
plane. Remind that for the massless case the pole exists at = — (2567r/A)^. In 
Fig. 19, we show the plot of H^ip^) for real p"^ (left) and the absolute value in the 
complex plane (right). 
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Fig. 19: Polarization operator for D=5 as a function of /m? for real J9^(left) 
and the absolute value in the complex plane (right). 



One can see that for negative the polarization operator is always greater than 1 
(in units of 3^^-^), for positive it is greater than 1/2 and then becomes complex. 
The absolute value in the complex plane is also always greater than 1/2. This means 
that depending on the value of dimensionless parameter ^ = A^m one has different 
possibilities: for ^ < 327r^ the pole exists at negative real p^, for 327r^ < ^ < 647r^ 
the pole exists at positive real < Am?. For ^ > 647r^ there are no poles at all. 
In this phase a theory is free from unphysical states. A similar situation, but in 4 



dimensions, was discussed in 



So, it looks like by choosing parameter ^ one can get rid of the unitarity problem. 
However, it reappears the other way. Indeed, one can see that the denominator of 
the propagator in this phase becomes negative. It is also negative at p^ = 0. In 
the scalar case the value of the a field propagator at = defines the effective 
potential of fields after integrating out the auxiliary field a. This way the negative 
value of the propagator leads to effective potential with negative quartic coupling, i.e. 
unbounded from below. In the case of the gauge theory the value of the denominator 
at = defines the sign of the residue of the gauge field propagator at = 0, i.e. 
the metric of the gauge field. Negative sign apparently leads to the "wrong" metric 
which is also not acceptable. Thus, the presence of a pole at the real axis is certainly 
a problem. 

The situation is different in D = 7 dimensions. Here one has 



4v/^(20 - 3a) - 3(a - Af ln( 



(7.6) 



128ya 

Notice the sign difference compared to the D = b case which means that here there 
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is no pole in the Euclidian region but in the complex plane. In Fig. 20, we present 
the same plots as above but for D=7. 



-4 -2 2 4 
_l I I I I I I I I I I I I I I \ I I I I I 




Fig. 20: Polarization operator for D=7 as a function of /m? for real p^(left) 
and the absolute value in the complex plane (right). 



One can see that the polarization operator in this case can take any value in- 
cluding negative ones. This means that complex conjugate poles exist for any value 
was already mentioned, they correspond to the ghost states and 
create trouble unless they are canceled. 

We now come to the last step of our analysis. According to Ref.[|l^, in the 
leading order the contribution of complex conjugated poles to a physical amplitude 
is canceled, thus preserving the unitarity in physical subspace. To check this, we 
consider the D=7 case and calculate the contribution from the conjugated ghost 
poles to the imaginary parts of the Feynman diagrams in the leading and next-to- 
leading order of the expansion. 

Consider first the one loop diagram shown in Fig. 21a. From the Kallen-Lehmann 
representation for the propagator of the auxiliary field ([7.3|) we take only the ghost 



terms ignoring the continuous spectrum. It corresponds to the following integral: 




a) b) 
Fig. 21: The leading order and next-to-leading order propagator diagrams 
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dk f R R* , 

+ .... ' (7-7) 



(k-p)^ \k^-M'^ k^-M'^\ 

where R and R* are complex numbers and M^, M^* are the masses of the conjugated 
ghost states. 

After integration, according to the dimensional regularization prescription, one 

gets 

R f dxil-xf'^p'x + M^+iTf'^ + R* f dx{l-xf/\p^x + M^-iTf/''. (7.8) 
Jo Jo 

For real p"^ the integrand apparently has no imaginary part being the sum of 
two complex conjugated expressions. The integration does not change this property: 
the contribution of the ghost states to the imaginary part (to physical amplitude) is 
canceled and only the continuous spectrum remains. What is crucial here is that the 
ghost states are conjugated having the opposite sign of the imaginary part and the 
same real part. 

Consider now the next-to-leading order diagram shown in Fig. 15b. In this dia- 
gram there are several ways how the ghosts might enter 

1. the inner propagator - non-ghosts, the outer propagator - ghosts; 

2. the inner propagator - ghosts, the outer propagator - non-ghosts; 

3. the inner propagator - ghosts, the outer propagator - ghosts. 

Consider the case when the ghost modes run in the inner propagator and in 
the outer propagator there is a continuous spectrum (non-ghosts). Then, using the 
integral representation (|7.3|) one has the expression 

A2(m2)3 + (81927r2)2 ^''^^ 
dkdq f R ^ R* 



(P - m?){{p - kfY{p -k-qf \q^ - q^ - M^* J ' 

Let us first take the integral over q. One has for the ghost part (and similar for 
the conjugated one) 

where we keep e to be finite since the two-loop integral diverges and omit the spectral 
integration over m^. The latter is real and is inessential. 
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R 

2e 



Evaluating the integral over k one gets 

R f\x{l-x))-^'^-'{l-xfdx C dy r dz{y - z)z-^/^+' 
Jo Jo Jo 

{ + T{2e){p^y{l - y) + m\l - y) + M^^—^z + iT^—^z)-^' 
- |/)V + 2/(m2 + ^r)(l - y)^x{M^ + tV)^] 
+ r(-l + 2e){p^y{l -y)+ m^l - y) + M^^^^z + ir^^^)^-^" 
[(7 -2e + 2)pV(l - yf + (7 - 2s)(M^ + ir)x] 

+ r(-2 + 2s){p^y{l -y)+ m\l - y) + M^^—^z + iT^—^zf-^' 

(7_2g)(3-2£) , 
^ 4 

Expanding over s one has singular and regular parts. The singular part is 

[\x{l - x))-^/\l - xfdx f dy r dz{y - z)z-^''' 
Jo Jo Jo 

{ + [{pY{1 - y)V + 2p\M^ + ir)(l - yfx{M^ + iVf] 
- -y) + m^(\ -y) + W^^—^z + iV^^—^z) [gpVfl - yf + liW + iV)x\ 

XX 

+ {p^y{l -y) + m^{l-y) + M^^—^z + iF^—^zfx^^}. 

X X S 

The remaining integrals over Feynman parameters are convergent and can be easily 
evaluated. One can see that due to the presence of iV the integrand is complex, but 
adding the complex conjugated term one gets the real polynomial of p^. 
As for the regular part, it contains 

log((p2y(l -y) + m^{l-y) + M^^—^z + iV^—^z)) 

X X 

and has a cut in momentum plane. However, the logarithm can always be presented 
in the form log(74e*''^) = log(A) + i(f), where both the modulus A and the phase 
depend on Feynman parameters. This means that adding the conjugated part one 
again gets the real integrand and, hence, the real function after integration. Here it 
is again crucial that the ghost states are conjugated and differ only by the sign of 
the imaginary part. 

Thus, wc conclude that the contribution from the conjugated ghost states to 
the imaginary part of the diagram is canceled and, therefore, the ghost states do not 
contribute to physical amplitudes. The same analysis can be carried out for the other 
choices of the ghost fields in Fig. 2 lb. Moreover, it seems to work in any diagram in 
all orders of the 1/N expansion since the reason for the cancellation is simple and 
obvious. This means that the unitarity in the physical sector is preserved. 
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The situation is somewhat similar to that in Ref . , where the mass generation 
problem was discussed in the context of higher derivative theory. Besides the physical 
states there exist non-physical states with a negative norm, but in the asymptotic 
states the negative norm excitations disappear thus preserving the unitarity of the 
theory. 



8. Conclusion 

We conclude that in higher dimensional scalar and gauge theories despite formal 
non-renormalizability it is possible to construct renormalizable expansion which 
obeys all the rules of a usual perturbation theory. The expansion parameter is dimen- 
sionless, the coupling is running logarithmically, all divergencies are absorbed into 
the renormalization of the wave function and the coupling. The original dimensionful 
coupling plays a role of a mass and is renormalized multiplicatively. Expansion over 
this coupling is singular and creates the usual nonrenormalizable terms. 

Properties of the 1/iV/ expansion do not depend on the space-time dimension 
if it is odd. In even dimensions our formulas after subtraction contain a logarithm 
which creates some technical problems in calculations but principally do not differ 
from the odd dimensions. 

Since the actual expansion parameter is dimensionless, all the Green functions 
get logarithmic radiative corrections and the cross-sections decrease with energy like 
in usual renormalizable theories without violating the unitarity limit. The running of 
the couplings depends on dimension and does not depend on Abelian or non-Abelian 
nature of a theory. This may be considered as a drawback of the expansion. 
Unfortunately, the preferable l/N^. expansion cannot be constructed in the same 
simple manner. 

We have demonstrated how one can deal with the problem of unitarity and 
unphysical pole states. The poles at the real axis can be removed by a proper choice 
of a dimensionless parameter ^ = XmPI"^'"^ which corresponds to the correct choice 
of the phase of a theory. However, this does not make a theory reliable. At the 
same time, the complex conjugated poles remain but fortunately their contribution 
to the physical amplitudes is canceled. We do not provide a rigorous proof of this 
cancellation but present the reason for it and several examples how it works in 
Feynman diagrams. Accepting this reasoning the theory seems to be unitary in 
physical subspace. 

We hope that this approach can be used in extra dimensional theories to get the 
scattering amplitudes. We expect that the behaviour of the cross-sections will differ 
from those of the Kaluza-Klein approach being closer to our approach based 
on the fixed points. 



Besides the already mentioned papers [0 there are several attempts to build 
renormalizable effective quantum gravity using a kind of expansion [03, where 
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the role of an expansion parameter is played by the number of space-time 
dimensions. The large D limit in this case is very similar to the large Nc planar 
diagram limit in the Yang-Mills theory considered by 't Hooft ||2^. The technique 
similar to the expansion is used also in [34|, where the author sums up the 

soft graviton corrections to the propagator of the scalar field and gets an improved 
propagator which decreases faster than any power of momenta. Though this partial 
resummation is similar to the expansion, the absence of an expansion parameter 
does not justify, to our mind, the selected set of diagrams. From this point of view the 
expansion is more consistent and contains the guiding line for such a selection. 
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